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Data from the first physics run at the Relativistic Heavy-Ion Collider at Brookhaven 
National Laboratory from the STAR experiment have been analyzed using three-pion 
\q . correlations to study whether pions are emitted independently at freezeout. We have made 
a high-statistics measurement of the three-pion correlation function and calculated the 
normalized three-particle correlator to obtain a quantitative measurement of the degree 
of chaoticity in the freeze-out environment. 
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Two-pion Hanbury Brown and Twiss (HBT) interferometry in principle provides a 
means of extracting the space-time evolution of the pion source at freeze-out produced 
in relativistic heavy-ion collisions |l|fj| . An underlying assumption of this method is that 
pions are produced from a completely chaotic source, i.e. a source in which the hadronized 
pions are created with random quantum particle production phases. Although two-pion 
HBT provides some estimation of this chaoticity, a better method is available using three- 
particle correlations. Normalizing the three-pion correlation function appropriately by the 
two-pion correlator, the effects from particle misidentification and decay contributions 
can be made to drop out thereby isolating possible coherence effects in the particle 
emission process. The resulting three-pion correlator r 3 provides the means of extracting 
the degree of source chaoticity by examining its value at zero relative momentum. 

2. Derivation of Chaoticity from Normalized Three-pion Correlator 

The measured observable in this analysis is the normalized three-pion correlator |3|]: 

r3 , Q ) = (C 3 (gg) ~ 1) ~ (C 2 (gig) - 1) - (C 2 (Q 23 ) - 1) - (C 2 (ga) - 1) (1) 
y/(C 2 (Q12) - 1) (C2 (Q23) - 1) (C2 (Q31) - 1) 



Here Q3 = \JQ\ 2 + Q23 + Q31 an d Qij = J — {pi— Pj) 2 are the standard invariant rel- 
ative momenta which can be computed for each pion triplet from the three measured 

momenta (p!,p 2 ,p 3 ). C 2 {pi,p 2 ) = H^iffe) and C 3 (Pi>P2,Ps) = ^^gf^^y , where 
P represents the momentum probability distribution. 



*For the full author list and acknowledgements, see Appendix "Collaborations" of this volume. 
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For fully chaotic sources r 3 approaches 2 as all relative momenta (and thus Q 3 ) go to 
zero. If the source is partially coherent, a relationship can be established between the 
limiting value of the three-pion correlator at Q3 = and the chaotic fraction e (0 < e < 1) 
of the single-particle spectrum ||: 

ir 3 (Q 3 =0) = v^^^. (2) 

Chaotic fraction e gives an upper limit on the value of the two-pion A parameter, which 
is sensitive to the number of coherent pairs in a sample. 

The three-boson correlation function C 3 (Q 3 ) is calculated from the data by taking the 
ratio anc ^ norma hzing it to unity at large Q 3 . Here A (Q 3 ) = is the three-pion 

distribution as a function of the invariant three-pion relative momentum, integrated over 
the total momentum of the pion triplet as well as all other relative momentum components. 
It is obtained by taking three pions from a single event, calculating Q 3 , and binning the 
results in a histogram. B (Q 3 ) is computed by taking a single pion from three separate 
events. 



3. Experimental Results and Discussion 

Data for the present results are from 1M events taken during the Year-1 physics run 
at STAR using the Time Projection Chamber (TPC) ||] as the primary tracking detec- 
tor. Two multiplicity classes were created by taking the 12% most central for the high 
multiplicity set, the next 20% most central for the low multiplicity set. For both multi- 
plicity bins, tracks were constrained to have pr in the range 0.125 < pr < 0.5 GeV/c, and 
pseudorapidity \rj\ < 1.0. In the range < Q 3 < 120MeV/c, approximately 150 million 
triplets were included in both the negative and positive pion studies. 

The C2 correlation function was corrected for Coulomb repulsion with a finite Gaussian 
source approximation The C 3 , correction factor is the product of three two-pion 
correction terms, obtained from the three pairs from the triplet. 

In calculating r 3 , the actual binned values of the correlation function for the various 
values of Q 3 are used instead of a fit ||. Q 3 , Q12, Q23 and Q31 are calculated from triplets 
of particles from the dataset and the three pairs that can be formed from the triplet. 
Eq. ([!]) is then evaluated (as a function of Q 3 ) using the binned two- and three-pion 
correlation functions, and averaged over the number of triplets in each Q 3 bins. 

The results for the three multiplicity bins are shown in Figures |1] (ir~) and |2| 
plotted as functions of Ql, and fitted to a function of the form: 

rs{Qs)=r 3 {0)-C 1 Ql-C 2 Qi. (3) 

This functional form is suggested by the theoretical analysis in || which shows that the 
leading relative momentum dependencies in the numerator and denominator of Eq. [I] are 
in even powers of Q 3 0. The fit range used is < Q 3 < 120MeV/c. 

The resulting intercepts ^(0) are shown in Fig. ^|, along with the results of WA98 and 
NA44. NA44 reported a result close to unity for Pb-Pb interactions, but a much lower 
result for S-Pb || , both with no clear Q 3 dependence. WA98 also reported a result close 
to unity at Q 3 = for Pb+Pb, and the (^-dependence in their result is similar to what 
we see in central collisions |§. 
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Figure 1. calculation for (a) central and 
(b) mid-central ir~ events. The fits shown 
use Eq. |3| to determine the intercept. Sta- 
tistical and statistical+systematic errors are 
shown. 



Figure 2. r% calculation for (a) central and 
mid-central ir + events. The fits shown use 
Eq. to determine the intercept. Statis- 
tical and statistical+systematic errors are 
shown. 



Figure |] shows the calculation of e for STAR's measurements, and for those from 
WA98 and NA44. The plot shows a systematic trend in the STAR results going from the 
peripheral bin to the central bin, with the central results showing a fully chaotic source. 
By calculating the number of coherent pair in the dataset from the chaotic fraction, a 
maximum value of the two-pion A factor has been determined which is consistent with 
the value of A found in the two-pion analysis HlO| . 

4. Conclusion 

In summary, we have presented three-pion HBT results for ^snn = 130 GeV data at 
STAR, and have shown that for the two multiplicity classes the STAR data indicate a large 
degree of chaoticity in the source at freeze-out. High statistics from STAR have allowed a 
normalized three-pion correlator calculation that extends to 120 MeV/c in Q 3 , and when 
used in conjunction with the formalism of partially coherent sources obtained from H, 
quantitative limits on the fraction of chaoticity are obtained which are in agreement with 
two-pion A parameter measured at STAR. STAR's measured values provide increased 
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Figure 3. Asymptotic value of r% from STAR 
and two other experiments ||[|. For STAR, 
central (a) and mid-central results are shown 
for 7r~ (circular markers) and tt + (square 
markers) data. For each STAR dataset, 
separate markers denote a fit using Eq. |3] 
(Quartic) and a fit to a quadratic equation 
(Quadratic). 



Figure 4. Chaotic fraction, calculated 
from Eq. |], and plotted for the same ex- 
periments as in Figure |3|. The data mark- 
ers are described in Figure |^. 



confidence in the validity of standard HBT analyses based on the assumption of a chaotic 
source. 
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